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ABSTRACT 


The reduction of circuit equations to normal form for numerical 
integration is considered for the general circuit where excess reactive 
elements, all types of linearly dependent sources, and nonlinear dissipa- 
tive and reactive elements are present. For the linear circuit, necessary 
and sufficient conditions for the existence of numerical solutions are 
considered and stated. For the nonlinear circuit, reduction to normal 
form is not always possible. Numerical solution is shown to be simpli- 
fied if certain a priori conditions are satisfied in formulating the original 
circuit equations. A new systematic reduction procedure is presented 
for obtaining the normal form equations. This procedure is also extended 
to a new procedure for obtaining transfer and immittance functions of the 


linear circuit from a proper tree formulation. 
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I, INTRODUCTION 

The formulation of circuit equations for their numerical solution is of 
current interest, The initial contribution was made by Bashkow, who 
formulated the concept of the proper tree, OO This was followed by 
the work of Bryant, and Kuh and Rohrer who refined the techniques. The 
object of their procedure, which was based primarily upon topology, 
is to obtain the circuit equations in matrix form so that one has a first- 
order differential equation. The approach uses the state-variable tech- 
nique, which has dominated the literature. However, it should be noted 
that other formulations of circuit equations, such as the TRAC Program, 
do not use state variables, 

Once the circuit equations have been formulated by the proper-tree 
concept, they have the following form 

Ax + Bw + Cz = Dx + Ew + FZ + Gu 
where x is the vector of voltages and currents associated with inde- 
pendent energy-storing elements, wis the vector of currents and voltages 
associated with dissipative elements, zis the vector of currents and 
voltages associated with excess reactive elements, uis the vector of 
independent sources, and the capital letters are coefficient matrices. 

For numerical solution it is necessary to alter the form to the normal 
a. x = Ax + Bu for the linear case and x = f(x,u,t) for the non- 
linear case, so that numerical integration can be applied directly. 

It has been shown that the circuit matrix equation, based on the 


formulation of the proper tree, gives the relation of the currents and 
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voltages of each circuit element to the currents and voltages of all other 
circuit elements. Of this set of circuit variables, currents and voltages 
of the circuit elements, there is a minimal set called states that are 
required to describe the circuit. If those variables which do not belong 
to the minimal set can be eliminated the resulting matrix equation will be 
of normal form, or an equation that can be readily converted to normal 
hOuM. 

Chapter II is a study of the circuit that contains linear circuit 
elements sree linearly dependent sources. The formation of the proper 
tree is discussed. A new systematic technique for matrix reduction 
that eliminates the unwanted circuit variables is developed. The circuit 
that contains no dependent sources has a matrix equation that always 
reduces to normal form and reveals that the number of states is equal 
to the number of independent energy-storing elements in the circuit. 

The introduction of dependent sources alters the coefficient matrices 

of the circuit equation. These alterations make the reduction process 
conditional and these conditions are discussed. In addition, it is found 
that the introduction of dependent sources may alter the number of states 
required to describe the circuit. Those reactive variables that are not 
destined to become states are designated as surplus. It is shown that 
the surplus reactive variables can always be eliminated from the 
equations, 

In Chapter [II the reduction techniques developed in Chapter II are 


utilized in establishing a new procedure for deriving driving-point 


impedance/admittance functions and transfer functions of the linear cir- 
cuit. This procedure enables a transfer immittance to be calculated 
directly from the proper tree formulation. Derivatives are replaced by 
their frequency-domain equivalents, 

In Chapter IV nonlinear dissipative elements are included with the 
linear circuit elements. It is shown that if normal form state equations 
are to be derived the currents and voltages of the nonlinear dissipative 
elements must be restricted in their relationship to the voltages and 
currents of mee circuit elements. This relationship may be controlled 
by a priori circuit modifications that increase the number of states. 

Chapter V extends the discussion to nonlinear reactive elements. 
Two new theorems relating to the conditions for the elimination of 
dissipative variables are stated. It is also revealed that the surplus 
reactive variables of the general nonlinear circuit must be associated 
linear elements and that these currents/voltages must not be related to 


the currents/voltages of nonlinear elements if normal form equations are 


desired, 


li, RESWeTehPOrcINEAR CIRCUIT EQUATIONS 
TO NORMAL FORM 


State equations are often used in circuit analysis. Normal form 
state equations, x = f(x,u,t), may be obtained from the matrix equation 


of the circuit's currents and voltages. 


A, FORMATION OF THE CIRCUIT MATRIX EQUATION 

In order to derive circuit equations in a consistent manner a proper 
tree is used. Any tree for the circuit may be used to obtain the circuit 
equations. However, if other than the proper tree is used, additional 
steps of transposing the variables to the proper side of the equation are 
necessary in the reduction process. The proper tree is drawn as follows. 
All voltage sources are designated as tree branches and all current 
sources are designated as links. This precludes a loop of voltage 
sources, which necessarily must be constrained so that the loop can be 
eliminated, or a cut-set of current sources, which must be constrained 
so that the cut-set can be eliminated, as illustrated in Figure 1. All 
capacitive elements, less the minimum required to prevent the formation 
of a loop, are designated as tree branches (C); those exluded become 
links (elastance, S). All inductive elements, less the minimum required 
to prevent the formation of a cut-set, are designated as links (L); those 
excluded become tree branches (inverse inductance, I). Resistive 
elements which do not complete loops are designated as tree branches 


(conductance, G); the excess become links (resistance, R). Thus the 
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order of precedence for drawing tree branches is (1) voltage sources, 
(2) capacitive elements, (3) dissipative elements, and (4) inductive 
elements. 

It is desirable to give consideration to elements whose currents 
and voltages are used as controlling terms for dependent sources. When 
a source is voltage dependent, the element across which the controlling 
voltage is developed should be a tree branch if this constraint does not 
upset the precedence for forming a proper tree. When a source is current 
controlled, ie element through which the controlling current flows should 
be a link. Should these arrangements not be possible, the controlline 
element's branch current may be written in terms of the link currents of 
elements forming a cut-set with the controlling element. Similarly, a 
link element's voltage may be written in terms of the branch voltages of 
elements forming a loop with the controlling element. Should the net- 
work be too complex to readily perform these substitutions and arrange- 
ments, this special consideration should be disregarded. Most circuits 
have very few dependent sources in comparison to other circuit elements, 
This consideration causes the dependency constants to appear in the 
circuit equations with the constants dictated by the network topology. 

Kirchhoff's voltage law is written for each fundamental loop defined 
by links L, R, and S. Kirchhoff's current law is written for each funda- 
mental cut-set defined by tree branches C, G, and T. All voltage sources 
appear in loop equations and all current sources appear in cut-set 


equations. Figures 2 and 3 illustrate the conceivable loops and cut-sets. 
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A loop defined by an L may contain dependent voltage sources, indepen- 
dent voltage sources, C's, G's andl's. A loop defined by an R may 
contain all of the foregoing with the exception of T's. A loop defined by 
an S does not contain G's andT's. Acut-set defined by a C may contain 
L's, R's, S's, and current sources. Cut-sets defined by G's do not 
contain S's and cut-sets defined by T's do not contain S's and R's. By 
using a proper tree other loops and cut-sets are prohibited. 

In the equations which follow, F (with appropriate subscripts 
relates link and tree branch voltages and link and tree branch voltages 
and link and tree branch currents according to the topology of the circuit. 


A prime designates the transpose of a matrix. and e. represent 


inde Ws S 


independent voltage sources. and Fe represent independent 


Ir ' Jo ' 
current sources. rand p (with appropriate subscripts) are submatrices of 
the proportionality constants of linearly dependent voltage sources. h 
and g (with appropriate subscripts) are the proportionality constants of 
the linearly dependent current sources. h and —~ are dimensionless where- 
as rand g have dimensions of ohms and mhos respectively. In the 

matrix equation which follows, fundamental loop voltage equations and 


fundamental cut-set current equations are written alternately as indicated. 


Thus, the circuit equations become (1). 
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B, THE REDUCTION PROCESS 

The reduction process includes the elimination of the variables 
associated with the dissipative elements, the transposition of the 
variables associated with the excess reactive elements, and the elimina- 
tion of the coefficient matrix from the left side of the reduced equations. 


The order of the first two steps is unimportant. 


13 


Partitioning (1) as indicated and lettering 





be] = vo i; [w= ; /- ai (2) 
€ | G hg 
enables (1) to be written as (3). 
wraaeg en = _ | 
iw 12 "13 | = ee tee Bie i. ie: 
: i 
a = a 
‘Aeuuee 2 ual Sr eh u w 
| ‘ | 
Zz 
Mel ec ae oo O32 223 | ne 


In (3), x designates the state variables, w designates the "unwanted" 
dissipative variables, and z designates the excess reactive element 
variables. Appendix A relates the submatrices of (3) to the submatrices 
oe) 
1, Elimination of Dissipative Variables 
The following rearrangement of (3) shows the relation of the 


variables, w, to the other variables and independent sources. 


_ E 
Pip 219) 271 213 | = Pay Py! re hea 

{ t \ ' ! 

} = — | wgei | | 

Pog 7 Po0) M1 = (221223) |, = Par Pos a 
! zi Es | oe 

P32 7 232) | 231 233] ES Bol Se iz. | 


_ — 


To have a determinate solution for the variables, w, the rank of 
the matrix multiplying w in (4) must be equal to the number of dissipative 
elements in the circuit. This matrix has dimensions n by m. Where n is 
the total number of circuit elements and m is the total number of 


dissipative elements in the circuit. Because of the dependent sources, 
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whose proportionality constants may be any real number, there exists the 
possibility that the matrix multiplying w in (4) will have rank less than 
m. In this event the solution for the variables, w, will be indeterminate. 

Let the rank of the matrix multiplying w in (4) be m-i because of 
row dependence. Rearranging the rows of (4) to group the rows such 


that the first tier is dependent upon the second tier leads to the following. 





aK ; ee le | 
py 12 a) ole Nyy >i ee |ecaihiimeni2 | a eis 
7 ‘ike C) , ! | 
M,, M,,.{w | Nv Nowe? oz PA ieee sk 4 e 
pale 22)| | (a1) enlnenc ee 22u eee 22 | aa 2 
(4a) 


M, > is a square matrix of rank m-i and dimensions m-i by m-i. The 
enclosed subscripts denote the dimensions of the w vector partitions. 
Since the first tier of M is a linear combination of the second tier, there 
must exist some matrix K s ha M j= [™ M a 

u i ome matrix K such t tk [M,, M,, 11 M12 


Multiplying the second tier of (4a) by -K and adding the results to the 


first tier yields 


Ir Bas -KN,,)- 1 
00 jiway | a KN,,) (N,,-KN,.) | x 
i | 
Moy My 01 n-i) [gen Moh: 
(T,-KT 51) @yy-KT,) | | : e Ke, | 
7 
TO Coy ee \ | “9 
L_ _ = . 


(4b) 
The second tier of (4b) yields the indeterminate solution for w. 
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Recalling that there were n elements in the circuit, m of which 


were dissipative, it can be seen that the first tier of (4b) contains n-m+ti 
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equations relating to n-m reactive variables. There is no guarantee that 
solutions obtained from a set of n-m of these equations would be con- 
sistent with solutions obtained from a different set. Assuming that the 
state equations obtained from the n-m+i equations were consistent, the 
substitution of these states and their derivatives into (4c) does not yield 
a determinate solution for the dissipative variables. Tnerefore, itis 
necessary that the coefficient matrix of the w variables in (4) have rank 
equal to the number of dissipative elements in the circuit if the w 
variables mic be eliminated. 

The submatrix ee ; has dimensionsm bym. If the rankism, 


ae 


then the solution for w is found in the second tier of (4). 


-] : . 
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Substituting (5) into the first and third tiers of (4) yields 





e ey a ' B | fx | V 

ce). (ELRTE. fi : 
OS alge u 5 || Me 

where 
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m = Ps) 4Q35-P39) (Py9-Qo) Pay (79) 
A= Pag *1Q307P a9) Ce as P03 (7h) 
= ox pe ioe ge a5 (73) 
Vo = ©, tQy9-P39) Poo Ooo) ew (73) 

In the event that (P,,-Q,,) is singular, the matrix inverse used 


in (5) would not exist. However, the n rows of (4) may be rearrange to 
provide a square matrix of rank m. This matrix may then be used in lieu 
227 92) in an equation similar to (5) yielding results similar to 


(6). Failure to achieve a suitable arrangement of the rows (4) leads to 


of (P..-Q 
an indeterminate solution for the variables, w. Since these variables 
are resistor currents and voltages and must remain finite, indeterminate 
solutions are not anticipated in circuit analysis problems. 
2. Transposition of Excess Reactive Variables 
To achieve the normal form equations from (6), z and Z must be 
transposed. This transposition yields 
x V 
| 
| : ~ ee (8) 


| 
ee 2 


== — —_— — 


— a 


pel a ¢ 
6 





m™ 


= 





The variables, Z, are currents and voltages of reactive elements 
that were deemed excess in the formation of the proper tree. Ina circuit 
free of dependent sources these variables would never become state 
variables. Because of the dependent sources (8) shows that this 


distinction can no longer be made. Thus, if 


vo = 





[" E 
(8) becomes 
Ay = By + U (10) 
Further reduction depends upon the rank of the matrix A. IfA 


is nonsingular, then (10) yields the normal form 


: -1 -l 
Vee Bt Dial) (11) 
If the matrix Ais singular, (11) is not valid and additional steps in the 


reduction process are required, 


3. Elimination of Surplus Reactive Variables 


The matrix Ain (10) can be singular because of linear dependence 
between the rows or between the columns or both. If the rows of (10) are 
regrouped forming one tier which consists of the independent rows of the 
matrix A and a second tier with the remaining rows and these tiers are 
partitioned to form square submatrices along the principal diagonals of 


A and B, the following results. 


Ay Ael [hl 4 
A A ’ j 
Mad Mid | Tb Pi Ra | 4 oe 

. | = ate (12) 
A A 
| di “dd | He Pai Pad | Gal U4 | 


In (12) the subscripts i and d denote independence and depen- 
dence respectively. 

Since the second tier of the matrix Ain (12) is a linear 
combination of the first tier, there exists some matrix K such that 


2 
A = | alee fees : | 
I an ag! KIA. Ea | . Multiplying the first tier by -K and adding the 
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results to the second tier yields 


Pe Aw) 3, | [oe B. 'y. lu. 
| hile Sate | i ic ii ke i a i (13) 
Ba | CaaS) ee aie: Yq “a 
Assuming that (Bg KB, is nonsingular, the second tier of (13) yields 
-] 
= = _ a = 
Yq (Bag7KBig) = (KB, Ba, ei a 3, 


Substituting (14) and the assoicated derivative equations into the first 


tier of (13) yields 


Axy. = SA CTO 1B Ue aes E*U, + Faun (15) 
where 
me + AB KB, 4) (SE 8 2a) (16a) 
ae - 6B Ct B, (Bag KB. |) (KB, -B a, (16b) 
oo B. , (By 47KB, ,) K a (16c) 
D* = -B, (8, 4-KB,,) (16d) 
EX = “A, (8 547KB, .)K (16e) 
ae = A, (8, ,-KB, 4) (1 6£) 


If A*¥ is singular because of row dependence the process used to 
derive (15) from (10) is repeated. Each repetition of this process may 
yield higher-order derivatives of the source terms. If A* is singular 
because of column dependence a process similar to that used to derive 
(15) is employed. 

Assume that the matrix Ain (10) is singular because of depen- 


dence between the columns. Regrouping and partitioning (10) 
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accordingly yield an equation identical to (12) where 





, = a | [K| (17) 
Multiplying the first tier of (12), partitioned according to column 


dependence, by as yields 







| | |, 

| 7 oie ee: (18) 
i :\ouune: ag we: B B | .y | U 

| di dd | bel di dd Gy _ d 


Multiplying the first tier of (18) by -A.. and adding the results to the 


di 


second tier yields (19). 





lt kK a wa: nae yd ADU 
aia : ipl ii id “8 ee 
} f=! ; + 
| -A._A -A_A taal -A_.A. 
ie - iXq Pai di ii Bi) Bad Gage el al Ua di ii Ys 
(19) 


Substituting the solution for Ya from the second tier of (19) 


and the associated derivative equations into the first tier yields 


KAS = * ee ok + * + * [ a eee 
AY BY: ae ra Eth Ud ee 
where 
See eee (Si avtp i (B_.-A a+ ) (Piha) 
dd “di ii “id di di ii ‘ii 
=i 7 Za zi =) 
* =_ _ -A pa 
: ij gee id'Pad a ae iQ (Bai di ae (2 1b) 
a = = Si 
*GC = A + A. -A A 21 
ii Bg Baa aii 3, a arte (21¢) 
ee 7 (pee A. fa) ) = (21d) 
id‘ dd “di id 
*E = =K{B. =A ays yea not (2le) 
dd “di ii “id di ii 
*F = K(B.,-A avtR Nee lt) 
dd “‘di ii ~id 
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If the matrix *A is singular and has row dependence, the 
process beginning with (12) is repeated. If there is only column depen- 
dence, the latter part of the process is repeated. 


In the foregoing it was assumed that the matrices (B ~KB, 4). 


dd 


=a 
ad “ALA, a B. qe used in (19) through 


used in (13) through (16), and (B 
(21), were nonsingular. If these matrices are singular, there will exist 
some rearrangement of the columns of (13) and (19) that will provide 


nonsingular matrices that can be used in place of (B “KB. ) and 


dd 


(B dd “AG AL, ii "Big 


This statement is proved by contradictions that show 
that the rank of the second-tier submatrix of y coefficients in (13) or 
(19) is equal to the number of variables to be eliminated. First, itis 
shown that there are no dependent rows; secondly, it is shown that the 
number of independent columns is equal to or greater than the number of 
variables to be eliminated. 

Let the following represent the arrangement of (13) or (19) in 
which the independent rows of the second tier of the y coefficient 


matrix have been separated from the dependent rows, the dependent rows 


being placed in a third tier. 











A 
ii id i Bai Pa oe cA 
| i == * * ; 1 | 4 * 
0 0 " Bs Bid | , Us (22) 
7d | 7d 
Le it It on ele it 
: : i i Pad pee 
a 1t tl = —_ uke ® 2 
Since |B Ai B id _| K [B*, Bd \. K being the matrix that shows 


the linear combination of the second tier of (22) that results in the third 


tier, multiplying the second tier by -K and adding the results to the third 


tier yields 





om A... | en, , we 4 
il id; -..- | tht rel _ vi 
A : ly, | 
= * 2 ae | * ®) 
1 0 0 ae Bai Bad | | US | (22a) 
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It can be seen that the third tier of (22a) contains only equations 
of constraint on independent sources. Since these sources cannot be 
constrained and remain independent, it must be assumed that the second 
tiers of the y coefficient matrices in (13) and (19) contain no dependent 
rows. 

Assume that the second tiers of (13) and (19) do not contain a 
sufficient number of independent columns in the y coefficient matrices 
to perform the required substitutions. Let the following represent the 
rearrangement of the second tiers of (13) or (19) which provides one 


square, nonSingular submatrix, the first tiers not being shown. 


Toe C6 ly | U*) 
, 4 y | | 
| Sa a = (22b) 
| O CA Cy5 4 : | U 
Let Ci be the nonsingular, square submatrix. Since the second set of 


columns of C is a linear combination of the first set, there must exist 


some matrix K such that 


wa petals 


C = Co, | [K | (226) 


Multiplying the first tier of (22b) by Ci 


yields 
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— 4 i aes. a ie 
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Multiplying the first tier of (22d) by “C,) and adding the results to the 


second tier yields 


lo ey is cy jut 
! 4 i222) 
0 0 yr Um Op Cent 


= 
a et Ue arang es 


oF 
Again constraint equations on the independent sources have appeared, 
Thus, it must be assumed that there are a sufficient number of inde- 
pendent columns in the second tiers of the y coefficient matrices in (13) 
and (19) to perform the operations required to eliminate the surplus 
reactive variables, Hence, the rank of the second-tier submatrix of y 


coefficients in (13) and (19) is equal to the number of variables to be 


eliminated, 


C. CONCLUSIONS 

In a circuit without dependent sources the number of states required 
to describe the circuit is equal to the number of capacitive branches plus 
the number of inductive links in the proper tree. When the circuit 
contains dependent sources, the proper tree does not necessarily 
Perecmine the number of states, It can be seen by comparing (7) to (1) 
that additional states maybe required if the dependent sources are 
controlled by currents or voltages of reactive elements made excess by 


the formation of the proper tree, 
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since the dependent source proportionality constants may take on 
any value, there exists the possibility that fewer states than predicted 
by the proper tree may be used to describe the circuit. Because of 
matrix singularities some of the reactive variables become surplus. The 
removal of these surplus variables introduces derivatives of the inde- 


pendent source terms and reduces the number of States. 


D, EXAMPLES 
iPeeeincredcee in States Duc tem eoendenmnclihees 
Consider the circuit of Figure 4 which illustrates the effect of 
dependent sources that requires additional states to describe the circuit. 


The circuit equations for Figure 4(a) are the following. 
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With the absence of dissipative elements (22) may be expressed in the 


form of (6) in which 
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(24) is in the form of (13) and the variables Vor and I. may 


be eliminated. With the dependent sources present as shown in Figure 


4(b). 
r BS C, 0 
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and the equations in the form of (24) become 
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If 6 is nonsingular, none of the variables can be eliminated. If 


6 is singular with rank of one, either Vo1 or IL-3 may be eliminated. 


Hence it can be seen that the presence of the dependent sources of 
Figure 4(b) caused an increase in the number of states required to 
Ges Chibe Lac eineule. 
Zee Deer a Sel DLatcamolle to Dependent Sounces 
Consider the circuit of Figure 5 which illustrates the effect of 
dependent sources that reduces the number of states by causing matrix 


Singularities, The circuit equations for Figure 5 are the following. 
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The equations of (27) are converted to the form of (3). Since 


th i = = = = 
ere are no excess reactive elements Pay Q3) Pao Q 
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From (28) it can be seen that (P -Q5 9) is singular and can 


Ze 
not be used in the solution for the dissipative variables. Noting that 


P Q, >) is nonsingular, the first two rows of (28) are used to solve for 


IQs 


these variables, 
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Substituting (29) into the second tier of (28) yields 
: ; , a = =n ie rz. : iv | (30) 
|Go1 G5CR, | Vo! [2 0 Vol ce 
(30) is in the form of (13). From the first row it can it can be 
seen that Vo = V. Thus substituting V for Vo in the second row of 


(30) yields the final result. 


Goll, = “L “G,CG,V (31) 


From (31) it can be seen that Vo is not needed to describe the 
circuit of Figure 5. Had the proportionality constant r been any other 
value, Vo would have been a state. Thus it can be seen that dependent 


sources may reduce the number of states required to describe the 


Cimeuit . 
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III, LINEAR CIRCUIT EQUATIONS 
REDUCED TO IMMITTANCE FORMS 


Linear circuit equations derived by the proper-tree concept can be 
reduced to impedance, admittance, and transfer functions if the circuit 
equations are augmented by output equations. When dealing with 
transfer/immittance functions the Laplace operator with zero initial con- 
ditions must be used in lieu of derivative notation. If the circuit equa- 
tions are x = Ax + Bu and the output equations are y = Cx + Du, then x 
may be found from the circuit equations and substituted into the output 
equations to yield Y(s) = [C(sI-A) B+D]U(s). This is an extended 
determination ofa ge immittance which reduces to the usual one 
when there is a Single input anda single output. This procedure requires 
one to obtain the state equations and output equations and from them 
derive equations that yield the desired functions. The procedure develop- 
ed here enables one to derive the desired functions directly from (1) by 


reduction techniques. 


A, DERIVATION OF DRIVING POINT ADMITTANCE/IM PEDANCE 

If a network is driven at one set of terminals by a voltage source 
and a proper tree is drawn of the network, only one fundamental cut-set 
will be described by the source. The sum of the currents of the link 
elements forming the cut-set with the source is the current passing 
through the source. Thus, the circuit equations may be augmented with 


an output equation which expresses the source current in terms of link 
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currents. Solving for the source current in terms of the source voltage 
results in an equation of the following form. 

nS) = ovo (32) 
I(s) is the source current, Y(s) is by definition the input admittance 
function, and V(s) is the source voltage function. If V(s) is unity (32) 
becomes 

I(s) = Y(s) (33) 

If the network is driven by a unit current source, only one funda- 

mental loop opi be defined by the source when the proper tree is drawn. 
The voltage developed by the source is the sum of the voltages of the 
branch elements forming the loop with the source. Thus, the circuit 
equations may be augmented with an output equation which expresses the 
source voltage in terms of the branch voltages. Solving for the source 
voltage in terms of the unit current source results in an equation of the 
following form. 

V(s) = G(s) (34) 
V(s) is the source voltage function, and Z(s) is the input impedance 
MUN C tHe mM" 


Writing the circuit equations in the form of (1), with LI replaced 


L 


by sLI., etc., and performing the substitutions given by (2) yields an 


L 


equation similar to (3), as follows. 
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Augmenting (35) with the output equation yields 
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where @ is the output current/voltage to be determined and Eg F ow! 
and ae are terms which relate the output to the circuit voltages and 
currents. The coupling terms, F* and F**, may include dependent 
source constants as well as topological factors. 

In the reduction process of Chapter II, to obtain the normal form 
equations, all dissipative currents and voltages were eliminated, The 
elimination process caused the dissipative parameters to appear in the 
coefficient matrices. The output equation identifies the currents and 
voltages of interest and to obtain an immittance the proper-tree equations 
should be reduced by eliminating those currents and voltages that do not 
appear in the output equation. Thus, after writing the equations in the 
form of (36), the unwanted variables are eliminated. The reduction 
process for these variables is identical to the process used to eliminate 
the dissipative variables in the derivation of normal form equations. 
(see Appendix C) 

The elimination of the unwanted variables, those for which the total 
output coupling (F** -F*) is zero, reduces the dimensions of the matrix 
equation. In addition, the reduction may cause the coefficients of the 
reduced equation to become polynomials or ratios of polynomials of the 


variables. The reduced equation is of the following form. 
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where the subscript r denotes the alteration of the coefficient matrices 
and y denotes the circuit variables that were not eliminated. Any of the 
elements of Pie Qe ay and aa may now be functions of the variable s. 
However, this notation has been omitted for simplicity in writing. 


Solving-for y in the first tier of (36a) yields 


st = fh-2,]7 [oy] oon 


Substituting this solution for y into the second tier of (36a) and solving 
for 6 yields 


7 ae 
= Kk = * aa 
9 [F2 al |P. Q, | vey (37) 
If 6 is the current of a unit voltage source, then the operations 


indicated by (37) yield an admittance function. If 6 is the voltage 


developed by a unit current source, (37) yields an impedance function. 


B, DERIVATION OF TRANSFER FUNCTIONS 

Network transfer functions may also be derived from the foregoing 
equations by letting @ become a voltage or current at terminals not 
directly connected to the source, Care must be exercized when the out- 
put is to be an open-circuit voltage as in the case of a transfer impedance 
function. There may be elements of the circuit that are in series with the 


Open-circuit terminals that do not affect the open-circuit voltage. These 


47) 


elements should be excluded from the proper tree and the resulting 
Cimemit equations. 

If the network is excited by a unit current source and @ is the open- 
circuit voltage at another set of terminals, then (37) yields a transfer 
impedance. If @is a current, the right side of (37) becomes a current 
gain function. If the network is excited by a unit voltage source and @ 
is a current through another set of terminals, the right side of (37) 
becomes a transfer admittance function. If @is a voltage, a voltage 
gain riawtise is derived. 

By letting 8 become an output vector, both input and transfer func- 
tions may be derived simultaneously. Thus (37) may be extended to 
multiple-port circuits with multiple excitations by identifying the sources. 
Since the output equations determine the unwanted variables, it may be 
necessary to perform the elimination process in stages to obtain various 


input and transfer functions. 


C. SUMMARY 
The steps to the derivation of driving-point and transfer functions 
are summariZed as follows: 
1, determine the sources needed to derive the desired functions 
2, eliminate the variables that appear in none of the output 
equations 
3. perform the operations indicated by (37) on those output equa- 
tions where the output is coupled to all the remaining variables 


including those variables that are not coupled but not effective. 
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D, EXAMPLES 
Consider the network shown in Figure 6a. Assume that the input 


admittance, Y for port l-1' is desired. Excite the circuit with a unit 


lal 


voltage source at port 1-1'. The proper tree is shown in Figure 6b. The 


augmented circuit equations are 
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The last equation is the output equation which shows the funda- 


mental cut-set relation of the source current and 1, The remaining 


il e 
circuit variables are to be eliminated. Elimination of the dissipative 


variables from (38) yields 
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Elimination of the capacitive variables from (39) yields 


sai] ] ] i = E 
sLI | (Ro oe) (J 1 
+. 
Jal | 2 sC, + Gy SC, + G. SC, G. Ei. 
i 1 
sLI | = ( ) = (—_.—__ ) | Lae 
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, =, a — as 7 a en | 
| 3 l 
| Yh i 0 0 
(40) 
Eliminating L» and performing the operations indicated by (37) 
yields 
° 3 Z 
+ + + 
Y _ oa a a5 d (41) 
Soe iS ee oS | encetet 
where 
a = L,C,°, (41a) 
= + 
b G,1,0, GAL,0, (41b) 
= + 
@) GGL, 7 (41c) 
d = Gy (41d) 
eS eS (41e) 
= + 
f G,L,L,C, +R,L,C,C, +G,1,L,C, (41£) 
g = G,G2L,L, + GR, LC, + R,G,L,0) + Lc, + Lc, 
aS, (41g) 
h = GR, G3L, + Git, + GL, + RoC, + GL, (41h) 
i = G.R,+1 (41i) 
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Assume that the two-port & parameter matrix for the network of 


Figure 6a is desired. 


current Sources is shown in Figure 6c. 
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From the output equations it can be seen that I 
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The proper tree with both ports being excited by 


The augmented circuit equations 
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Should be eliminated, The result of their elimination and the relocation 
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Performing the steps necessary to obtain the form of (37) yields 





ih 
sC.+G 
1 i 
= SL, 
1 ljfl sL| 0 
_ L,C,s +L,G, +] 
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i l 
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Completion of the multiplication indicated by (44) reveals the desired 


parameters. ,The results are 





gue jaa ls (45) 


2 
= = a ee LC (46) 
3 12 g21 422 2 l 
L,C,s -- L,G,s + ] 
The solutions have been left in a form that allows one to check their 
accuracy by inspection. 


The hybrid parameters could be derived by exciting the circuit at 


one port with a voltage source and a current source at the other port. 


IV. REDUCTION TO NORMAL FORM 
WITH NONLINEAR DISSIPATORS 


Because of the nonlinear current-to-voltage relations of semi- 
conductors, the numerical solutions for circuit containing these devices 
may become complex if the state variable approach is used. The circuit 
can be represented by a set of normal-form equations that are augmented 
by equations which also must be solved to obtain the nonlinear currents 
and voltages that appear in the normal-form equations. The general 


form of these equations are 


es ee Cf, W) - DE, (i) 
v= A,x + Bou + f,(i,v) (47) 
i= Ax + Bu + f,i,v) 


where the f'x are nonlinear functions and v and i are the voltages and 
currents of the semiconductors. If f, and fy are zero, then (47) reduces 
to one nonlinear differential matrix equation which may be integrated 
directly. 

Diodes either by themselves, or as part of the eee 


equivalent for the transistor may be considered in terms of their current/ 


voltage characteristics. For the solid-state diode it is possible to 


write 
oa kv 
i, = I (e ad -)) 
and 
i. +I 
teakes © d O 
ee lat ; ) 


The exact nonlinear relation is unimportant for the present consideration 
so that the form 


i = f(v,) (48a) 


and the inverse form 


Vv, = g(i,) (48b) 


are used here. 

Werther and pao and Selah” consider the diode to be a 
current-dependent voltage source or a voltage-dependent current source 
where the controlling voltage or current is written in terms of the voltages 
and currents of other circuit elements as in (47). This treatment requires 
the solution of network equations for the determination of the controlling 
voltages and currents. If the both the current function, (48a), and the 
voltage function, (48b), are single valued (the tunnel diode is exluded), 
then the nonlinear dissipative element may be treated as the linear 
dissipative elements are treated when writing circuit equations. That 
is, the device may be used as either a branch element or a link element 


when forming a proper tree, 


A, TOPOLOGICAL CONSIDERATIONS 

The nonlinear elements are given precedence between capacitive 
elements and linear dissipative elements when establishing the branches 
of the proper tree. When writing a voltage equation for a loop defined by 
a nonlinear, dissipative link element, a term on the left side of the equa- 


Hom such as RI, is replaced by the voltage v The term corresponding to 


qd’ 
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I, of the right side of the equation is f(v ,) . Similarly, when writing a 
current equation of a cut=set defined by a nonlinear, dissipative element, 
the term on the left side of the equation is Ly and the term on the right side 


is g(i,). The matrix equation for the circuit will be similar to that of (1). 


B, ELIMINATION OF DISSIPATIVE VARIABLES 

If state equations are to be obtained the dissipative variables must 
be eliminated. While the linear dissipative variables may readily be 
eliminated, the elimination of the nonlinear variables depends upon the 
relationship of the currents/voltages of the nonlinear elements to the 
voltages and currents of other circuit elements. This may be seen by 
the following discussion. 

1. Nonlinear Variables As General Functions of Other Variables 

Assume that the nonlinear voltages/currents are functions of all 

other voltages and currents. After writing the circuit equations and 
completing the substitutions given by (2) the equations become 


So. 


fu Pia Pag Pra Pas] 4 Pa @12 913 S14 ‘i ot es. 
Par Poa Pog Pag Pas . wa Qo 223 Sag P25) | * ee 
gi (32 “33 3a 35 w — |@sy “a9 33 * 44 a5 WaT ley 
Par Paz Pas Pag Pas’ | Ya yy Say P43 gq Sas] |E%g) “¢ 
P51 ?52 P53 Psa Pos | ta! es, 252 253 254 S55 [Fa] Ja) 
(49) 
where x,zZ, and w are given by (2) and P and Qi, (i,j, = 1,2,3) are 
given by Appendix A, Py and Qi G=4,5;3,=1,2,3,4,5 and j = 4,5; 
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i=1,2,3,4,5) are matrices which relate the voltages of nonlinear 
dissipative links to other voltages in the fundamental loops which are 
defined by nonlinear elements, and the currents of nonlinear branches to 
the currents of other elements which are in the cut-sets that are defined 
by nonlinear, dissipative branch elements. 

The first two tiers and corresponding columns of (49) may be 
combined by transposing z and z. Thus, by forming y from x and z as 


was done in (9), (49) becomes the following. 


y a fe 
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As was discussed in Chapter II, the w variables may be eliminat- 
ed if the rank of the combined coefficient matrices for these variables is 
equal to the number of these variables, Their elimination reduces (50) 


to the following. 


Sq Sip ue Y Ma ig Be ee " 
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The multiplication of (51) by [K] yields (47). Further 
reduction depends upon the elements of Mo: M,.1 M,. , and M,.. This 


may be seen by writing the second two tiers of (51) as follows. 


43 S24 36 onl eben 3 M39 M33 4 
Equation (53) suggest several undesirable solutions for the 
nonlinear currents and voltages, Transcendental solutions may be 
required. If+these are not required, solutions which reveal that a non- 
linear current/voltage is a nonlinear function of other currents/voltages 
may result. The substitution of the latter solutions into the first tier 


of (51) yields functions of functions. If M,,.=M,,=M,,=M,, = 0, 


these undesirable solutions cannot exist. By comparing (52) with (50), 


it can be seen that the desired conditions prevail if By 3B, .=B,5=B, 2 Bay 


= =B,,=B,, = 0. This implies that the fundamental loops defined 

42 43 44 
by nonlinear elements must not contain other dissipative elements, 
either linear or nonlinear. Similarly, cut-sets defined by nonlinear 
elements must not contain other dissipative elements. These restric- 
tions must also apply to dependent sources which are controlled by the 
currents/voltages of dissipative elements. Thus, if desirable equations 


are to result, the nonlinear, dissipative variables must be restricted to 


being functions of reactive variables only. 


2. Nonlinear Variables As Functions of Reactive Variables 
If the currents and voltages of nonlinear elements are not allowed 
to be controlling terms for dependent sources and these currents/voltages 
can only be related to capacitor voltages or inductor currents, then (50) 


becomes the following. 
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The eliminations of the w variables yields the following. 
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The elimination of V4 and i, from (55) yields 
Ky = M14 + ue + M,,{(M,.¥ + V,) + M,,9(M, ae) (56) 


If Ki in (56) is nonsingular, then (56) readily yields the normal 


form equations. If Rid is singular some of the variables in y must be 
eliminated. The elimination of the surplus reactive variables may be 
possible but very difficult, Assume that the left side of (56) is zero 
and that one of the variables y is to be written in terms of the other 
variables and sources, It can be seen that the solution of a transcen- 
dental equation may be involved. To investigate the conditions neces- 


sary for reduction when K iS Singularaeciveuit thatwsoemee Or 


iol 


dependent sources is considered, 


For a circuit that does not contain dependent sources the circuit 


equations reduce to the following. 
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The matrices Qa 4s Q.. ae) and Qe will be other than zero 


Anes 
whenever a fundamental cut-set is formed linear dissipative elements 
and nonlinear dissipative elements. Likewise, Qys and Qe will be 
other than zero whenever a fundamental cut-set is formed with nonlinear 
dissipative elements. The elimination of the linear dissipative variables 


) and g(i .) in the reduced equations 


gives rise to coefficients for fv, 4 


causing the solutions for Vy and 1g to become transcendental. If 1g is 


eliminated from (57), then Qas and Qe4 cause the solution for V4 to 
become transcendental, Thus, normal-form state equations will not 
result from the reduction of circuit equations when fundamental cut-sets 
are formed with nonlinear dissipative elements and other dissipative 
elements. 

If the matrices mentioned in the preceeding paragraph were Zero, 
the nonlinear voltages/currents would be functions of the state variables 
and sources only and all dissipative variables could be eliminated. The 


resulting normal form equations would be as follows. 


x= Ax + BU+CU + Df(Ex + FU) + Gg(Hx + JU) (58) 
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C. CONCLUSIONS 

The circuit equations for circuits containing nonlinear dissipative 
elements do not readily reduce to normal-form state equations. Bya 
priori circuit modifications of placing an inductor in series with those 
nonlinear elements which have loop equations formed with variables that 
are not destined to become state variables and placing a capacitor in 
parallel with those nonlinear elements which have cut-set equations 
formed with variables that are not destined to become state variables, 
the nonlinear elements' voltages and currents are forced to become state 
variables. While these modifications do give rise to additional states 


they allow the circuit equations to be reduced to normal-form equations. 


D. EXAMPLE 
Consider the circuit shown in Figure 7a. The circuit equations are 


the following. 


CV “cL OG) Vo | 0 0 

LI ee ee : 1 0 v.| 
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This equation cannot be put in normal form because of the coupling 
between diodes. If the circuit is modified as in Figure 7b, the circuit 


equations become 


45 


a ™ | | 

nel | 0 Clie | | Vou 0 ~ 7 

| ies i at : 1, + | l o - 

| Vat ee ee Egy) | 0 " ae) 
a 1 0 Lo 0 | [Ey go) | ol wiley 


The diodes are now uncoupled and (60) may be placed in normal 


ieleaal 
CV. “G, 0 1. Vo IO) _ 
° | pW o)) | | 
= . - x ' 
CVG] Terme a! | + 
— foci 
LI PI Ve Shs 0 Oo; Aa 
— _ | LL ib | L 
(61) 


(61) is obtained by eliminating v , from (60). Normal-form equations 


d 
are obtained by multiplying (61) by 


=| 


| 


G 
© 
© 


(62) 


lo ‘@) 
—) @) 

— 
jo (ex 
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V. REDUCTION TO NORMAL FORM 
WITH GENERAL NONLINEAR ELEMENTS 


When nonlinear inductors and capacitors are present, the circuit 


equations may be reduced to the following form. 


| x | £(x) | o- 
| | xe x | 
i 
| 
a lwo > w \ + e ; (63) 
i) a ij | , a 
| | ! | 
volves loi): Piro 
. a meee) [od | 
where 
: oo O.) 
x 4 i £,@) =ly | =e): ( ) (64a ,b) 
ol cree 
Ve fs@3 
; £, (2) =|, ° E (_) (65a,b) 
Tr 
B. is the flux associated with inductors and a. is the charge 


associated with capacitors, p, and Q. are functionally related to the 
inductor currents and capacitor voltages as given by (64b). For the 
linear case g ~ LL and Oe (65a,b) give the general relationships 
for excess reactive elements. 


If x and z in (49) are replaced by f(x) and f(z) respectively, the 


result is (63). Since z and f(z) may be transposed, it is possible to 


form 
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£ (x) | 


y= and f(y) =|j~ |. (66a,b) 
rae 


se | 


Ne 


If fy) is used in place of y in the equations that begin with (50), the 
arguments for reduction with nonlinear inductors and capacitors would 
be obtained. Thus, if the dissipative variables are to be eliminated, 
the conditions for the elimination of these variables from (49) must also 
apply to (63). 

The conditions for the elimination of the dissipative variables are 


¢ 


stated in the following theorems. 


Theorem 1. Ifa coefficient submatrix for the linear dissipative variables 
in the equations which relate these variables to the other circuit 
variables and sources has rank equal to the number of linear dissipative 
elements, the linear dissipative variables may be eliminated and normal- 


form equations may be derived. (For proof see Chapter IIB1) 


Theorem 2, If each of the voltages or currents of nonlinear dissipative 
elements, nonlinear dissipative variables, are linearly related to only 
independent sources and the reactive voltages and currents that are 
functions of variables that are destined to become state variables, then 
the nonlinear dissipative variables may be eliminated and normal form 
equations may be derived. (For proof see Chapter IVB2) 

Having eliminated the dissipative variables from (63), the resulting 
equations are as follows, 


Ay = Bi ty) th 1) at ee pew) fame) #9 Ee eee, + Uc 
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If Ain (67) is singular there exist surplus states to be eliminated. 


Since these states do not appear in functions f, and om which result 


d 
from the elimination of the nonlinear dissipative variables, an equation 


of the following form must be solved to eliminate the surplus state 


functions fy) and the surplus state derivatives. 


fice —™ 
i ty), | 
E se 
0 | M N £ Wg i) (68) 


— 


where the subscripts i and s denote the state functions to be retained 
and the surplus state functions respectively. While the surplus state 
functions EY). may be written in terms of the other state functions and 
sources, the solution for the surplus state derivatives involves partial 


derivatives of both the surplus functions and the functions to be retained. 


df (y) , df (y), , 
mY s_ y = N ly ee y, + N ly (69) 
oY. S 2A il 


From (69) it can be seen that if the surplus states are to be elimin- 
ated the surplus state functions must be linear and must be related to 
retained functions that are also linear. Otherwise, the partial deriva- 
tives must be entered into the coefficient matrix of the retained 
derivatives and the resulting equations become nonlinear differential 


equations that cannot be placed in normal form, 


APPENDIX A 


General Submatrix Relations to Circuit Equation Submatrices 
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APPENDIX B 





pe 37 V3" ¥4 
(a) Voltage Source Loops 





(b) Current Source Cut—Sets 


FIG ELIMINATING VOLTAGE SOURCE LOOPS 
ANDO CURRENT SUURCE CUT—SiErs 
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FIG 6. EXAMPLES 3 & 4. IMMITTANCE FORMS 
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(b) Modified Circuit 


FIG 7 EXAMPLE 5. INTRODUCTION OF ADDITIONAL 
STATE TO OBTAIN NORMAL FORM REDUCTION 
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APPENDIX © 


Elimination of Unwanted Variables 


Consider the following equation and the elimination of the variable 


A A (ae F.9) B B B x e | 


- _ | 
7 6 ae Le ee we 





ile 32> celia aoe ce Baa] Py 122 | 


If (A__-B...) is nonsingular, the second equation may be written as 


me (+8 4 eB F(z) -A 


21 He He ae 


++ 
Da” oo) 


Za 
Substituting this value for w into the first and third equations and 


combining terms yields 


. 7 . =e 
lil 
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ime pattern ior the multiplication and Combining of terms is readily 


established. 
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